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Abstract 



The tunnehng between two parahel two-dimensional electron gases has 
been investigated as a function of temperature T, carrier density n, and the 
applied perpendicular magnetic field B. In zero magnetic field the equilibrium 
resonant lineshape is Lorentzian, refiecting the Lorentzian form of the spectral 
functions within each layer. From the width of the tunneling resonance the 
lifetime of the electrons within a 2DEG has been measured as a function of 
n and T, giving information about the density dependence of the electron- 
impurity scattering and the temperature dependence of the electron-electron 
scattering. In a magnetic field there is a general suppression of equilibrium 
tunneling for fields above B = 0.6 T. A gap in the tunneling density of states 
has been measured over a wide range of magnetic fields and filling factors, 
and various theoretical predictions have been examined. In a strong magnetic 
field, when there is only one partially filled Landau level in each layer, the 
temperature dependence of the conductance characteristics has been modeled 
with a double- Gaussian spectral density. 
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I. INTRODUCTION 



In addition to conserving energy E, resonant tunneling between two parallel two- 
dimensional electron gases (2DEGs) required the conservation of the in-plane momentum 
hk. In zero magnetic field, the conservation of both quantities allows the investigation of the 
broadening of the electronic states within an individual 2DEG, in a way that is not possible 
with conventional transport measurements. For example, recent experimental studies! have 
shown that the lifetime of the two-dimensional electrons can be measured directly by tunnel- 
ing spectroscopy. The lifetime broadening is characterized by the electron spectral function 
A{k, E), which describes the probability that an electron in a particular fc-state has energy 
E. In general, A{k, E) is strongly peaked near the free-particle energy Ek = h^k'^/2m*, with 
a width r determined by the lifetime r of the electron. 

In a strong perpendicular magnetic field, when the Landau level filling factor u = hn/eB 
is less than unity, it is assumed that the electron traverses the tunnel barrier between the 
two 2DEGs in a much shorter timescale than the charge rearrangement around the injected 
electron (and around the hole left behind in the emitter layer). The energies of the injection 
and extraction processes create a gap in the tunneling density of states (DOS). Many-body 
theories predictETO that for a fixed filling factor, u < 1, the width of the energy gap is 
approximately 0.4e^/47reZs (where Ib = \JU/eB is the magnetic length), or if the electron 
liquid forms a regular latticei the gap energy is e^/47rea (where a = l/y/rvK is the inter- 
electron spacing within a layer) . At fixed filling factor both predicted forms for the gap have 
the same a/B dependence with magnetic field. In an earlier paper we reported! tunneling 
characteristics in the high magnetic field regime {u < 1). The theoretical predictions do 
not match the measured linear field dependence, A ^ OAAhuJc (where huc = heB/m* is 
the cyclotron energy), that we have observed!'! in high magnetic fields both at fixed filling 
factor and constant carrier density. 

From tunneling measurements of 2D-2D samples with different barrier widths, Eisenstein 
et alM argue that the high field gap energy is largely determined by the in-plane Coulomb 
energy, modified by an interlayer excitonic binding energy. The form of the gap used to fit 
their data was = /(i^) e^/47rea — OAe"^ /Aned, where /(z/) is some universal function of 
filling factor, and d is the quantum well center-to-center separation. Although the interpre- 
tation of our data is different to that of Eisenstein et alB (there is no universal /(z/) which 
allows us to plot our high field data in this form), we measure gaps of comparable energy in 
lower mobility, but otherwise similar GaAs/AlGaAs double-layer structures. 

In a 2D-3D tunneling system Ashoori et a/.0 measured a filling factor independent 
suppression of tunneling for > 1. The temperature dependence of the conductance was 
fitted to a model with a linear variation of the DOS about the Fermi energy Ep. The 
half-width at half-depth of the gap was measured to be linear in magnetic field, with a 
dependence Aa ~ huJc/20. Taking into account the difference in the two tunneling systems 
and the definitions of Aa and A, the two gap widths can be arguedS to be comparable. 

Here we consolidate earlier workpiili at zero and high magnetic fields, and present data 
taken at intermediate fields. The rest of this paper is organized as follows. The samples and 
measurement details are described in Sec. The basic tunneling formalism is reviewed in 
Sec. |ITT| , showing how the tunneling conductance is determined by the underlying electron 
spectral functions within each 2DEG. In Sec. |V|we present tunneling measurements in zero 
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magnetic field, investigating the various scattering mechanisms which are responsible for the 
linewidth of the resonance. The tunneling properties in a perpendicular magnetic field are 
investigated in Sees. |V] and |VT[ 

II. EXPERIMENTAL DETAILS 

The measurements presented in this paper were obtained from two samples, C751 and 
C770, each consisting of two modulation-doped 180 A wide GaAs quantum wells separated 
by a 100 A (C770) or 125 A (C751) Alo.33Gao.67As barrier. The double-well structures 
were grown by molecular beam epitaxy on an n"*" epilayer, which had previously been pat- 
terned by a focussed Ga ion beam to form backgates.EJ Optical lithography was used to 
fabricate a mesa, Ohmic contacts, and to define surface Schottky gates. Small front (back) 
gates near the Ohmic contacts were used to selectively deplete the upper (lower) 2DEG, 
forming an independent contact to the lower (upper) layer. The two 2DEGs overlap in a 
100 /imxlSO /im tunneling region. The as-grown carrier densities of the upper and lower 
layers in the tunneling region are ni ~ 3 x 10^^ cm~^ and n2 ~ 2 x 10^^ cm~^. The corre- 
sponding, as-grown, low-temperature, mobilities of the two layers are /ii ~ 8 x 10^ cm^/Vs 
and /i2 ~ 2 X 10^ cm^/Vs. The carrier densities in the tunneling region are independently 
controlled by voltages Vg-^ and Vg^ applied to large area front- and backgates. 

Two-terminal conductance measurements of the sequential depopulation of Landau lev- 
els (LLs) in each layer were used to determine the relationship between V^^ ^ ^''^d the 
corresponding carrier density ni.2. For example, in C770 the top layer has a density 

= {r.lVgjV + 2.9) X 10^^ cm-2, where -0.4 < Vg, < 0.35 V. The carrier density at 

= can vary by as much as 10% from one sample cooldown to another. By contrast, the 
constant of proportionality (7.1) between Vg^ and rii is always the same, and is in good agree- 
ment with a model where the top 2DEG and the front gate form a parallel plate capacitor 
of separation 1000 A. 

The differential tunneling conductance G = dl/dVsd was measured using constant voltage 
(0.1 mV) ac lockin techniques as a function of carrier density n, temperature T, and applied 
dc interlayer source-drain voltage Vgd- For measurements in a perpendicular magnetic field 
B the carrier densities were always matched, and the given values of u refer to the filling 
factors in each layer. 

III. TUNNELING FORMALISM 

Tunneling accesses the single-particle DOS (the DOS available for the addition of an- 
other electron), compared to conventional transport measurements which provide informa- 
tion about the DOS available for the excitation of an electron from the Fermi sea into 
a conducting state. Zheng and MacDonald" showed theoretically that 2D-2D tunneling 
characteristics in zero magnetic field can be used to determine the spectral function of an 
electron. Within layer i [i = 1,2 (upper, lower)] scattering introduces an energy broadening 
Ti = h/2Ti, where is the scattering time. In the Born approximation the spectral function 
within each 2DEG will have the Lorentzian formS 
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where E is measured from the bottom of the 2D subband. Using standard notation^! the 
tunnehng current /(V^d) between the two layers is calculated from the two spectral functions 
to be 

J(\/,d) oc J (fk^ J d'k2 |Tfc,,fc, |2 J dE 
A,{ki,E-Eo,)A2{k2,E-Eo,) 

x[/(E-/xi,T)-/(E-/i2,T)], (2) 

where f{E,T) = 1/(1 + exp{E/kBT)) is the Fermi distribution function. The chemical 
potentials /ij and the 2D subband energies Eq- are defined in Fig. |I|. 

The integral is simplified by the following assumptions: the absolute conservation of 
momentum (such that the interlayer tunneling matrix element T^^ is proportional to 
(^fci.fcz)) low temperature (/cb^ ^ T, with Fermi functions approximated by step functions), 
and weak disorder (F ^ Ep). Using these assumptions, substitution of Lorentzian spectral 
functions into Eq. gives the tunneling current 

G F^ 

where AE-p = E-p^ — Ep^^, and Go is a constant with the dimensions of conductance. The 
tunneling linewidth is the sum, F = F1+F2, of the widths of the Lorentzian spectral functions 
in the two layers." 

In this paper we report differential conductance G = dl/dVsd measurements. The zero- 
bias conductance obtained by differentiating Eq. |^ is 

G(Vsd = 0)(x Jd^kAi{k,EFjA2{k,Ep,). (4) 

Using Lorentzian spectral functions, the differential conductance as one subband edge (con- 
trolled by a gate voltage V^) is swept past the other i 

Equation |^ predicts Lorentzian G{Vg) characteristics with a linewidth F and peaked with 
magnitude Gq at AEp = 0, when the carrier densities in the two wells are equal. 

We are interested in the differential conductance as a function of V^d, taken when the 
carrier densities have been matched at V^d = 0. However, upon application of an interlayer dc 
voltage the capacitance between the two layers unavoidably mismatches the carrier densities, 
and it is incorrect to set AEp = in the derivative of Eq. ^ to obtain the conductance 

There are in fact two contributions to AEp: One component AEp is due to the applied gate 
voltages (this is zero when the carrier densities are matched at V^d = 0), and the other, AEp , 
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is a capacitive effect. AEp is determined by the change of carrier density Sn = ±CKd/e 
induced by V^d, where C = e/d is the capacitance per unit area between the two 2DEGs. 
The capacitive change in AEp is AEp = —2 \ 5n\ /po = —PeVgd, where po = rn* /irh?' is 
the 2D density of states, and (3 = a-B/2d (where the Bohr radius is ob ~ 100 A in GaAs). 
Substitution of AEp = AE^ + AE^ = AE^ - /JeKd into Eq. |, followed by differentiation, 
gives the modified lineshape 

Setting AEp = (for matched carrier densities at V^d = 0) gives the lineshape 

G(K.)^G„I!(I1^(£M1^. (8) 

(r2 + (ei4i(i - im 

which has a linewidth of r/(l — /?), in contrast to the linewidth F given in Eq. ^ A 
typical value of (3 is 0.15. Therefore, linewidths obtained from G{Vsd) measurements are ap- 
proximately 18% larger than those obtained from equilibrium conductance measurements.ll3 
Similar capacitive corrections should also be applied to linewidths obtainedl from plots of 
I/Vsd versus l^sd- 

We can only quantify capacitive effects at i? = 0, though the magnitude of such effects 
in a perpendicular magnetic field can be estimated by replacing po with ps, the density of 
states at the Fermi level in the outer Landau level. If the outer LL is half-filled, Pb ^ Po, 
then the capacitive change of the Fermi energies AEp will be negligible. The opposite 
result is obtained near integral filling factors and so we have not made any conclusions from 
tunneling measurements in this regime. 



In Sec. we model high field tunneling measurements obtained at filling factor 
z/ = 1/2. To model the data we have guessed the form of the spectral density Ai{E) = 
l/(27r)^/ d'^k Ai{k, E), and have calculated the tunneling current using the expression 

I{Vsd) <xJdEAi{E- EoJ A2{E - EqJ 

x[/(E-pi,r)-/(E-p2,T)]. (9) 



IV. TUNNELING IN ZERO MAGNETIC FIELD 

Figure ^ shows the zero-field differential tunneling conductance G as a function of V^^ , 
the front gate voltage; by fixing Vg^ , the carrier density in the bottom layer was held constant 
at n2 = 3.25 x 10^^ cm~^. The equilibrium measurements were taken at T = 3.0 K (squares) 
and 19.0 K (circles). Due to a misalignment between the frontgate and backgate in the 
tunneling region, there is an unwanted second resonance (which is independent of Vg^), as 
well as a small constant background conductance. In Fig. |]both of these weak contributions 
have been removed, and the corrected data at the two temperatures have been fitted to the 
Lorentzian lineshape 

G{Vg,) = ^ 2, (10) 

" 1+ ^ " ^ 



5 



which is Eq. prewritten exphcitly in terms of the front gate voltage V^i. The fits, shown as 
sohd hnes, are with Vq = 0.046 V, SV^ = 0.0194 V, and Gq = 16.1 /xS at 3 K, and with Vq = 
0.045 V, SVg = 0.0360 V, and Gq = 9.9 /iS at 19 K. We have also tried fitting the resonances 
to a Gaussian and to the derivative of a Fermi function (/' ~ sech^[(Vgj — Vo)/(5V^]). Of 
the three functional forms, a Gaussian fit is the poorest, /' is considerably better, but the 
overall shape and the tails of the resonance are best fit by a Lorentzian. The half width 
at half maximum of the low-temperature Lorentzian resonance, SVg ~ 19 mV, corresponds 
to an energy broadening F ^ 0.5 meV. This linewidth F is the sum of the widths of the 
spectral functions in the double-layer system, F = h/2Ti + h/2T2- We relate the linewidth 
to the scattering rate, defining the average scattering rate of the double-layer system to be 

Non-equilibrium measurements provide an alternative method for determining the 
linewidth F. Figure ^ shows a comparison of the non-equilibrium G(Kd) (squares) and 
equilibrium G(V^J characteristics (circles, and offset vertically) for the same 77,2 for sample 
C751. To demonstrate the validity of Eq. |^, the non-equilibrium data was obtained with a 
small mismatch of carrier densities, causing a slight shift of the central G{Vsd) peak to the 
left of Vsd = 0. The solid line through the G'(Kd) data is a fit to Eq. with AE^ = 0.08 meV, 
corresponding to a carrier density mismatch of 0.7%. Though small, this mismatch is enough 
to cause a noticeable asymmetry in the G{Vsd) characteristics. The fits to both the equi- 
librium and non-equilibrium traces in Fig. ^ give F = 0.384 meV and Gq = 3.03 fiS. The 
fitted value /3 = 0.14 ± 0.01 is obtained from the non-equilibrium trace and corresponds to 
an interlayer separation of d = 365 =b 30 A, somewhat larger than the quantum well center- 
to-center separation of 305 A deduced from the growth parameters. However, both parallel 
magnetic field measurements (using the method of Ref. |l^ and self-consistent calculations 
suggest a center-to-center wavefunction separation of d ~ 340-350 A. 

The successful fits to both the equilibrium and non-equilibrium resonances in Fig. |^ 
provide strong evidence for a Lorentzian spectral function in zero magnetic field. We have 
been able to fit Lorentzian lineshapes to the equilibrium tunneling resonances obtained over 
a wide range of carrier densities, 0.3 < n < 3.2 x 10^^ cm~^, at temperatures from 3 K to 
19 K. From the n and T dependence of the linewidth F, we have investigated the electron- 
impurity, electron-phonon, and electron-electron contributions to the overall scattering rate 

r/n = T-' = + r-p\ + T-l (11) 

Figure ^ shows the temperature dependence of the equilibrium linewidth F at matched 
carrier densities of n = 0.91, 1.62, 2.19, 3.04 x 10^^ cm~^. The solid lines are fits of the form 
F = a+bT"^; at high n the fits are good, but at lower densities the measured linewidth departs 
from quadratic behavior. The fitting parameters are both functions of n, with typical values 
a ^ 0.5-0.8 meV and 6 ~ 5 x 10^^ meV/K^. Only the electron-phonon and electron-electron 
scattering terms in Eq. ^ are expected to show significant temperature dependence. In the 
temperature range considered here, the sheet resistivity pxx of each 2DEG shows only a weak 
temperature dependence, indicating that the mobility is dominated by impurity scattering, 
and that, by comparison, the temperature dependent electron-phonon contribution is small. 
We therefore associate the temperature dependent component of the tunneling linewidth 
with that due to electron-electron scattering Fe_e. The electron-electron scattering rate in a 
clean 2D metal is given by the well-known theoretical expressionEl 
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h Ef f kBT\ ^ + In P'^TF I _^ ^ 



(12) 



where gxF = 2/aB is the 2D Thomas- Fermi screening wavevector. At fixed carrier density 
Ti/tc-c exhibits an approximate temperature dependence, in agreement with the sohd fine 
fits used in Fig. ||. 

At low temperatures, the electron-electron scattering is negligible, and is dominated 
by electron-impurity scattering. Figure ^ shows a log-log plot of the equilibrium linewidth 
versus carrier density at T = 3 K. The fitted gradient is a; = —0.47 ± 0.05, which is 
approximately the power-law behavior expectedli^ for small-angle scattering from remote 
ionized impurities. 



^c-imp 2,771*3 V 2 



^n-'/^ (13) 



where no is the sheet density of impurities, and s is the spacer layer thickness. Equation |13 
is valid for 4kps ^ 1, where kp is the Fermi wavevector. In our samples the ionized donors 
are distributed over distances of 400-700 A from the center of each quantum well. Taking the 
smallest distance, s = 400 A, the condition Akps ^ 1 is easily satisfied for n = 3 x 10^^ cm~^, 
though the applicability of Eq. |1^ becomes questionable at the lowest carrier densities. 

From experiments over a range of carrier densities and temperatures we have fitted over 
300 measured equilibrium linewidths r(n, T) to the expression 

r(n, T) = re_in,p + Fe-e = n""'^ + X [Eq. 12]. (14) 

If F is measured in units of meV, n in units of 10^^ cm~^, and T in K, the best fit was 
obtained with the numerical values ai = 0.77 ± 0.01 and a2 = 3.06 ± 0.09. The dashed lines 
in Fig. ^ are obtained from this best fit. 

The magnitude of the electron-impi|irity scattering rate, Fe_imp = ain~^'^ in Eq. |l^, is 



— — n " t-.-iiiij-' . ^ — |— -1 7 — 

^ ^..^ , il but is fifteen times smaller than the theoretical 

prediction of Eq. 0. Discrepancies between theory and experiment have previously been 
resolved^ by including correlations between scatterers. In the high n, low T regime, the 
ratio t^'^/t^^ ~ 7 (where Tt is the mobility lifetime) agrees with measurements^ of similar 
mobility single- layer GaAs/AlGaAs heterostructures. Furthermore, the scattering time r ~ 
1.2 ps is comparable to the Dingle time 1.1 ps (obtained from Shubnikov-de Haas oscillations 
measured on a Hall bar made from the same wafer), confirming the importance of small-angle 
scattering in determining the tunneling linewidth. 

To look at just the electron-electron scattering term. Fig. |] shows Fe-eZ-EFJ'lotted versus 
k-QT/ Ep. This representation of the data is expected to be nearly universal,a though there 
is some dependence on carrier density inside the square bracket of Eq. [1^. This weak 
dependence on n is illustrated by the two solid line traces of 3.06x[Eq. |T2[ (our best fit for 
Fe-e) for n = 0.3 and 3 x 10^^ cm~^ (upper and lower solid curves, respectively). Equation ^ 
is plotted for n = 1.3 x 10^^ cm~^ and without any prefactor (a2 = 1) as the dashed line 
labelled GQ in Fig. |^. Since the original derivation by Giuliani and Quinnlil there has 
been some confusion about the value of a2 that multiplies this expression. Fukuyama and 
AbrahamJni calculated the expression 
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(15) 

and this is shown in Fig. ^ as the curve labelled FA.@ More recent theorieJ^lii highlight 
the mathematical errors in Ref. |18| and suggest that due to the nature of the tunneling 
experiment, there are both electron and hole contributions to the linewidth. These two 
considerations result in 02 = 7r^/2 = 4.9, giving the third dashed line in Fig. ^ labelled 
JM. Without any energy dependent corrections the various predictions (GQ, FA, and JM) 
can be compared to our equilibrium measurement a2 = 3.06. The difference between the 
predictions of GQ and JM and our fit is purely in the value of the prefactor 02; whereas the 
curve FA has a functional form which does not describe the data very well. 

The greatest discrepancies between theory and experiment in Fig. ^ occur for large 
k^T/Ep, for which there are various possible reasons. First, it is assumed in the derivationOJ 
of Eq. |12| that Ep ^ k^T, a condition which breaks down at low densities and high temper- 
atures. Second, when the matched carrier densities are reduced there is an increase in the 
disorder of the two layers, and the 2DEGs may no longer be considered as clean 2D metals. 
The increased level of disorder may also render the Born approximation invalid, in which 
case the spectral functions and the tunneling resonance are no longer Lorentzian. Third, the 
derivation of Lorentzian G(V^J characteristics (Eq. |iy) from Lorentzian spectral functions 
is only valid when F -C Ep; the two energies become comparable when n = 0.5 x 10^^ cm~^. 
The fourth reason is experimental: the pinch-off characteristics of the gate (V^J distort the 
symmetric equilibrium lineshape, increasing the uncertainty of measured linewidths at low 
n. 

In our analysis of the temperature dependence of the equilibrium linewidth we have 
followed an approach similar to that of Murphy et al. ,1 where the non-equilibrium linewidths 
were measured from I /Vsd versus V^d characteristics of double 2DEG structures of higher 
mobility (/i ~ 10^ cm^/Vs). In our samples the measured electron-electron scattering rate 
is enhanced over Eq. |12| by a factor 02 ~ 3, which is two times smaller than that measured 
in Ref. 0. The two measurements are different, and non-equilibrium linewidths are expected 
to be larger than equilibrium linewidths for two reasons: an 18% enhancement can be 
attributed to charge transfer effects (see Sec. ^), and a further 15% (Ref. ^5|) to 40% 
(Ref. ^6]) enhancement may be introduced by the energy dependence of Fe-e. In order to 
achieve an enhancement factor 02 ~ 6 close to that of Murphy et a/.i recent theorieJ^lii 
have concentrated on the energy dependence of Fe_e, but have neglected the interlayer charge 
transfer effects. In Fig. |^ we show that, in our samples, we can obtain good agreement 
between equilibrium and non-equilibrium linewidths by including only charge transfer effects. 
Moreover, in further investigations over a range of matched and mismatched carrier densities, 
we find that the energy dependent corrections are consistently small. 

There are two major differences between our samples and those of Ref. ^ that might 
account for the different measured values of 02. The first is the level of disorder; in our 
samples the electron-impurity contribution to the linewidth is Fe-imp ~ 0.5 meV compared 
to 0.25 meV in Ref. 0. This greater disorder probably does not affect the applicability 
of Eq. 12, but will make the energy dependent corrections to Fe_e less important. § The 
second difference is the thickness of the tunnel barrier, which will affect the strength of the 
interlayer electron-electron interactions. For samples with quantum well widths of 200 A 
and barrier thicknesses in the range 175-340 A ( givmg d = 375-540 A), Murphy et a/.i 
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measure an enhancement a2 = 6.3 which appears to be independent of barrier thickness. 
In our samples the interlayer electron-electron distance {d ~ 340 A) is comparable to the 
lowest used in Ref. 0. This could suggest that it is the greater disorder in our samples that 
gives rise to the smaller enhancement for Fe-e. By way of comparison, we point out that 
Berk et al^ have measured electron-electron scattering rates similar to our own (in the 
units used here, 02 = 3), in 2DEGs with a similarly low mobility fi ~ 50-400 x 10^ cm^/Vs 
but more strongly coupled {d = 180 A). The effect of stronger tunneling (as quantified by 
the symmetric-antisymmetric gap Asas) on Fe-e has been calculated by Slutzky et a/.ii 



V. TUNNELING IN A PERPENDICULAR MAGNETIC FIELD 

The application of a perpendicular magnetic field B has a marked effect on the resonant 
tunneling properties between two 2DEGs. Figure |^ shows the equilibrium conductance 
G{Vsd = 0), at a matched carrier density of n = 1.0 x 10^^ cm~^, as a function of magnetic 
field at 0.1 K (solid line) and 1.3 K (dashed line). As B is increased, the oscillations 
in the tunneling conductance mirror the formation of LLs in both 2DEGs. However, in 
contrast to SdH oscillations, which are an in-plane transport effect, the oscillations shown 
in Fig. originate from the current flowing perpendicular to the 2DEGs.i2l The oscillations 
in G{Vsd = 0) can be interpreted in terms of variations of pb, the DOS at the Fermi level in 
a field B. We will describe the temperature dependence of the equilibrium tunneling in two 
magnetic field regimes, B < 0.6 T and B > 0.6 T. 

Below 0.6 T, the formation of LLs modulates the DOS, creating minima (maxima) at 
integral (half-integral) filling factors. The equilibrium tunneling amplitude depends on the 
DOS in both layers and shows quantum oscillations commensurate with pB- When the 
temperature is increased from 0.1 to 1.3 K, the low field tunneling conductance follows the 
behavior of p^; the amplitude of the conductance oscillations decrease as the modulation in 
Pb becomes weaker, and G(V^d = 0) tends towards its zero- field value. 

For B < 0.6 T the magnitude of the conductance maxima increase with magnetic field, 
following the expected increase in the degeneracy of a half-filled LL. However, for B > 0.6 T 
the amplitude of the tunneling maxima decrease with increasing magnetic field; this can be 
interpreted as a magnetic-field-induced suppression of the equilibrium tunneling between the 
two 2DEGs, with the formation of a gap in the tunneling DOS. When the temperature is 
raised to 1.3 K the suppression of the conductance is reduced, indicating thermal activation 
across the gap. The Fig. |^ inset shows low-field equilibrium tunneling for a matched carrier 
density of n = 3 x 10^^ cm~^. At this higher carrier density the crossover from SdH-like to 
activated behavior again occurs at -B = 0.6 T. These measurements at different temperatures 
show that the signature of the magnetic-field-induced gap is evident at filling factors as high 
as 1/ ~ 18. Clearly, above 0.6 T a mechanism other than the simple formation of LLs 
is needed to explain the unusual temperature dependence observed at half-integral filling 
factors. 

The formation of a magnetic-field-induced gap can also be investigated by non- 
equilibrium measurements, where at low magnetic fields a dip around zero bias is observed 
in the GiVsd) characteristics. Figure | shows G{Vsd) characteristics at various magnetic 
fields for a matched carrier density of n = 0.95 x 10^^ cm~^. In the trace at i? = 0.4 T 
there is a slight dip at Vsd = 0; with increasing magnetic field both the depth of the dip and 
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the separation of the surrounding peaks become larger. To characterize the tunnehng gap, 
we define A to be the voltage separation of the conductance peaks either side of V^d = 0. 
Previously, we have measure H^sd) characteristics for z/ < 1 and have successfully fitted 
our low voltage data to the expression I = Jq exp(— A/el^d)- At high fields the values of A 
obtained from such fits are precisely the same as those measured from the peak separations. i 
Although at low magnetic fields we cannot fit our data to this theoretical expression for /, 
we nevertheless use the peak separation as a straightforward measure of the gap. To test 
various theoretical predictions we have investigated A at fixed n and fixed z/, both as a 
function magnetic field B. 

Figure ^ shows the gap parameter A, measured from traces similar to those in Fig. |^, 
over a wider range of B. Due to the presence of edge states we disregard measurements close 
to the integer quantum Hall (QH) regime, as well as those near the fractional QH state at 
1^ = 2/3. The linear fit through the remaining points (solid circles) is A = (0.45±0.02)/ico'c — 
(0.19±0.05) meV. Figure ^shows that the gap is most cleanly observed at half-integral filling 
factors, and in measurements at fixed n there are more data points at low u than at high u. 
With this proviso the linear fit in Fig. ^ spans the whole field regime, though it is difficult 
to identify deviations at low B. We have also performed conductance measurements at a 
higher matched carrier density of = 1.48 x 10^^ cm~^, and Fig. shows the evolution of 
the the low field tunneling characteristics at this carrier density presented as a grey-scale 
plot. The figure is made up of G{Vsd) characteristics, similar to those in Fig. |^, taken at 
101 magnetic fields from S = to 2 T. To accentuate the data each trace has been scaled 
between (black) and 1 (white). The light structures around the center, marked with solid 
lines, are the two tunneling peaks surrounding the gap, with a splitting that is linear in 
magnetic field, A = (0.43 ± 0.03)^u;c — (0.26 ± 0.03) meV. Within experimental error the 
gaps measured at the two carrier densities, 0.95 x 10^^ cm^^ and 1.48 x 10^^ cm~^, have the 
same magnetic field dependence. 

To obtain more information about the gap at higher filling factors. Fig. |TT| shows mea- 
surements of Aj, taken as a function of B while maintaining a half-integral filling factor in the 
tunneling region. Gap measurements were obtained at z/ = 1/2, 3/2, . . . , 11/2 for seventeen 
matched carrier densities in the range 0.36 < n < 3.10 x 10^^ cm~^. At a given magnetic 
field, the higher the half-integral filling factor, the smaller the gap. The DOS is not clearly 
spin-split below B ^ 1 T, and so data points for the highest filling factor and lowest fields 
are subject to some error. The straight line 0.45^ti;c — 0.19 meV is included in Fig. 11 to 
demonstrate the linear B dependence of the tunneling gap at u = 1/2, and the sublinear 
nature of the gap at higher filling factors. At low carrier densities n < 1.5 x 10^^ cm~^, 
measurements of Ai,>3/2 converge to a value similar to that measured at z/ = 1/2. For the 
higher carrier concentrations the gap for u > 3/2 is reduced in comparison with A1/2. This 
observation, plus the scarcity of data points at higher filling factors in measurements at 
fixed n accounts for the same measured gap at 0.95 x 10^^ cm~^ and 1.48 x 10^^ cm~^. We 
stress that, as previously reported,0 the high field gap shows the same linear behavior both 
at fixed n and fixed u < 1. 

Figures ^ and |10| also show tunneling between LLs of different index. At high magnetic 
fields, the first inter-LL tunneling transition has been observedBS as a peak in the /(V^d) 
characteristics at eVsd ~ 1.3^ci;c, and the equivalent conductance measurement]! exhibits a 
peak at eV^d ~ 1.16^c- Both measurements show inter-LL transitions occuring at higher 
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than expected bias voltages, and it is proposedEI that this enhanced LL spacing is due to 
an effective mass reduced by many-body effectsU The arrows in Fig. |^ and the dashed hues 
superimposed on Fig. |T0| show the position of the first inter-LL tunnehng peaks either side 
of zero bias. Over the range 0.5 < 5 < 3 T the conductance peaks occur at eV^d = (1-07 ± 
0.03)^u;c, an enhancement that is only slightly less than that observed at higher magnetic 
fields. If the enhanced LL spacing is due to a reduced effective mass, our measurements 
suggest that the mass is reduced at higher filling factors by an amount much greater than 
that predicted by Smith et a/H 



Comparison with Theory 

There has been much interest in tunneling in a weak magnetic field where, even for filling 
factors z/ ^ 1, it is believed that there is an energy cost associated with the injection and 
extraction of electrons into and out of a 2DEG. It has been theoretically shownil that at low 
magnetic fields the bare Coulomb interaction (which typically exceeds the cyclotron energy 
huc) is screened by the filled LLs, and the renormalized interaction is sufficiently weak that 
it does not mix LLs. Therefore only interactions between electrons in the outer partially 
filled LL need to be considered; this is the starting point for many of the theories in the 
weak field regime. 

Following the treatment at high magnetic fields, the gap at lower fields and higher filling 
factors can be considered to be the Coulomb energy, where the cyclotron length Ic = Tik-p/eB 
replaces the magnetic length Ib as the relevant length scale. The semiclassical energy gap is 

where n is the interaction parameter, and is given by = a/a^ = \/2/a-B,k-p. Levitov and 
ShytovEil consider the problem of an electron orbiting in a ring of thickness Ib and radius Ic 
in one layer, tunneling into a similar ring in the second layer, leaving behind a hole in the 
first ring. The electrostatic energy of this two-ring capacitor creates an energy gap that is 
proportional to Eq. jTB], but multiplied by a logarithmic screening term. Alternatively, using 
both a hydrodynamic approach!! and a Hartree-Fock calculation,!! Aleiner and coworkers 
have obtained a gap of the form 

A2 = ^lnM). (17) 

It has been shownS that the total gap has contributions from two terms similar in form to 
Ai and A2, where Ai is the dominant term. 

In a recent theory, Fogler et a/.S propose that under certain conditions the ground state 
of the 2DEG in a weak magnetic field is a charge density wave (CDW) superimposed on a 
uniform background of filled LLs. Tunneling into a partially filled LL disturbs some of the 
correlations of the CDW, and the energy cost of this disruption gives rise to a tunneling gap 

"■^^""^ In f 1 + + A2 ^ 0.07^^e + A2. (18) 



y27r V r, 
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As B is decreased there is a predictedUl crossover from the spin-resolved QH regime (spin- 
split LLs) where the ground state is a CDW, to the spin-unresolved QH regime (well sepa- 
rated, but spin-degenerate LLs). The CDW ground state can only exist in the former regime, 
while at lower fields the ground state is a uniform electron liquid which is expected@ to have 
a gap of a different origin. In our samples the SdH oscillations of the individual 2DEGs show 
that the spin-unresolved QH regime occurs over a narrow range of magnetic field. The Fig. |^ 
inset shows that at = 3 x 10^^ cm~^ it is possible to observe a suppression of tunneling 
(above 0.6 T) in the absence of spin-splitting (which becomes distinct for B > 1.5 T). More- 
over, in all our measurements over a wide range of n, we see no clear change in the character 
of the gap once snin-splitting has set in. The B dependence of the low field gaps predicted 
by Aleiner et a/.@ and Fogler et a/.ii are determined by A2, which varies as A2 ~ 52 at 
fixed carrier density, and as A2 ~ -B at fixed filling factor. Figure p| shows that at fixed n 



the gap is linear (not quadratic) in B, and Fig. |TT] shows that at fixed u the weak field gap 
^u>3/2 is sublinear in B. 

The semiclassical gap (Eq. |1^) behaves as Ai ~ i? at fixed carrier concentration, and 
as Ai ~ \^ at fixed filling factor. For fixed n the measured gaps are indeed linear in B, 
though we find that the best fits are obtained with a negative intercept. Similarly the best 
fits of A,y to a field dependence are also obtained with a negative intercept, using the 



expression = a' yB/u — b'. In the inset to Fig. |Tl] we have plotted (Aj^ + y')v^l'^ versus 

V > 5/2. Taking b' = 0.4 meV the data collapse onto a single straight line with 
slope a' = 2.0. We conclude from measurements at constant n and constant u that the high 
filling factor gap is best described by 



e2 



A = (0.46 ± 0.03)^^ - (0.4 ± 0.1) meV. (19) 



The negative intercept could be interpreted as evidence for an excitonic effect (similar to 
that observed! at u = 1/2), and with a magnitude which is approximately independent of 
filling factor for u > 5/2. We were unable to obtain good agreement between Eq. ^ and the 
measured gaps at z/ = 1/2 and 3/2, instead much larger negative intercepts {b' = 4.2 and 
1.7 meV, respectively) were required. With our samples we have obtained! better agreement 
at u = 1/2 using a linear rather square root magnetic field dependence for the gap A1/2 up 
to 25 T. 



VI. TEMPERATURE STUDIES IN HIGH FIELDS 

In this section we use Eq. ^ to model high field G{Vsd) measurements taken at z/ = 
1/2, results which have been presented elsewhere.0 Figure |12Ka) shows the conductance 
characteristics at -B = 8 T at a matched carrier density of rii = n2 = 0.97 x 10^^ cm~^. 
As the temperature is increased from T = 1.5 K to 6 K, the equilibrium conductance 
G{Vsd = 0) increases, while the height of the surrounding peaks is reduced. Figure |l2|(b) 
shows G{Vsd = 0) in an Arrhenius plot, revealing activated behavior at low temperatures 
with an activation energy of E^, = 0.35 ± 0.01 meV. Above the activation temperature 
(pa 4 K) the equilibrium conductance departs from exponential behavior, and for T > 10 K 
the conductance decreases. We have measured at different matched carrier densities 
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and magnetic fields, while maintaining a filling factor of z/ = 1/2 in the tunneling region. 



Figure |13| shows Eg^ as a function of magnetic field, with a fit which shows that the activation 
energy is proportional to B. At u = 1/2 the ratio E^,/ A is 0.047±0.005, and is approximately 
independent of B. Activation plots of the high field gap were first obtained by Eisenstein 



et al, and for comparison we measure from Fig. 3 of Ref. |31] the ratio E^,/ A = 0.07 at 
5 = 13 T. 

Temperature can be used to probe the shape of the tunneling DOS, and previously 
Ashoori et a/.0 have used a DOS with a linear dip at the Fermi level to model 2D-3D 
tunneling measurements. A linear DOS does not give activated behavior, and to obtain 
such characteristics we have modeled the data using the double- Gaussian spectral density 

where the Fermi level Ep lies in the minimum between the two peaks. This tunneling DOS 
is a modified version of Eq. 14 in Ref. ^ which was originally proposed for z/ 3> 1 and 
Eq ^ k^T. By comparison with the original formula, we have introduced an adjustable 
parameter a in the denominator of the exponential. As T ^ the double-Gaussian spectral 
density consists of two 5-functions at £" = ±i?o, and so the gap A — > AEq at T = 0. At low 
temperatures the activation energy Eg, — > E^/a, and so A and E^^ share the same magnetic 
field dependence (that of Eq), in agreement with experiment. 



To simulate the conductance measurements in Fig. |12] we have chosen a = 5.3 so that 
the theoretical ratio E^/A = l/4a agrees with the experimental value of 0.047. /(V^d) 
characteristics have been calculated using Eq. |2y in Eq. ^, and by numerical differentiation 
of these characteristics we obtain the G{Vsd) traces shown in Fig. |l^(a). The value of 
Eq = 2.6 meV best reproduces the conductance and activation characteristics (see Fig. [T^(b)) 
at 8 T. From similar simulations of u = 1/2 data at different magnetic fields we have 



measured the B dependence of Eq; the results are displayed in Fig. |T^, with a fit showing 
Eq with a linear dependence on the magnetic field, Eq = 0.2huJc- 

The calculated G{Vsd) curves reproduce many features of the experimental data. First, 
there are bias voltages which delineate regions of positive and negative dG/dT; second, 
the conductance maxima move closer together with increasing temperature; third, it shows 
that at higher temperatures G{Vsd = 0) is no longer activated, but reduces with increasing 
temperature. There are, however, some differences between the model and experiment. The 
very low temperature limit of the experimental characteristics differ from those of the model; 
in the experimental system disorder introduces a finite width to the DOS at T = 0. We have 
tried to simulate this disorder by introducing a finite offset temperature, that is by replacing 
T with T + To in Eq. ^ where Tq is a constant. However, there was little improvement in 
the shape of the simulated characteristics. 



VII. CONCLUSIONS 

We have used both the equilibrium and non-equilibrium tunneling characteristics to 
show that the underlying spectral function, A{k, E), within each of the 2DEGs of a double- 
layer system is a Lorentzian. Equilibrium measurements show that the Lorentzian tunnel- 
ing linewidth has a carrier density and temperature dependence that can be attributed to 
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electron-impurity and electron- electron scattering. The measured electron-impurity scatter- 
ing rate is similar in magnitude to the Dingle time, and comparison with theory shows that 
correlations between scatterers need to be taken into account to obtain better agreement 
with experiment. The mathematical form, but not the magnitude, of the electron-electron 
scattering rate measured from the equilibrium lineshape is well described by theory. 

In a perpendicular magnetic field the equilibrium tunneling conductance oscillates as a 
function of B, reflecting the formation of LLs. For B > 0.6 T the magnetic field creates 
a gap in the tunneling DOS which suppresses the equilibrium tunneling. For fixed filling 
factor V — 1/2 the tunneling gap has a linear magnetic field dependence, whereas at higher 
half-integral fiUing factors the gap is sublinear and is best described by a semiclassical gap. 
We have modeled the temperature dependence of the conductance characteristics in a strong 
magnetic field using a double-Gaussian spectral density. 
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FIGURES 

FIG. 1. Conduction band profile of tlie double 2DEG structure, defining the various energies 



used in Sec. Ill 



FIG. 2. Equilibrium tunneling conductance G{Vsd = 0) as a function of the gate voltage 
V^^ controlling the carrier density in the upper layer, when the lower layer density is fixed at 
n2 = 3.25 X 10^^ cm~^. The traces were taken at T = 3 K (squares) and 19 K (circles), and the 



solid lines are the fits to Eq. 10 



FIG. 3. Comparison of non-equilibrium (squares) and equilibrium (circles) tunneling lineshapes 
at 0.1 K. The solid lines are fits to Eqs. ^ and|l^, respectively. For clarity not all the experimental 
data points are shown, and the equilibrium data and fit have been vertically offset by 1 ^S. 

FIG. 4. Temperature dependence of the equilibrium linewidth F 
for m = 722 = 0.91, 1.62, 2.19, 3.04 x 10^^ cm'^. The solid lines are fits to the form a + bT'^, 



and the dashed lines are obtained from the best fit to Eq. 14. 



FIG. 5. Log-log plot of the linewidth F versus carrier density n at T = 3 K. 

FIG. 6. Tn^/Ep plotted versus kBT/Ep. The dashed curves are the theoretical predictions 
of GqJiI FA,ii and JmII plotted for n = 1.3 x 10^^ cm"^. The sohd curves are the best fit to 
Eq. [l^ with a prefactor of 02 = 3.06, and with n = 0.3 and 3 x 10^^ cm~^ (upper and lower curves, 
respectively) . 

FIG. 7. Main figure: equilibrium tunneling conductance at a matched carrier density of 
n = 1.0 X 10^^ cm~^, taken as a function of perpendicular magnetic field. The solid line shows 
data taken at T = 0.1 K, and the dashed line shows data taken at T = 1.3 K. Inset: similar 
measurements taken at n = 3.0 x 10^^ cm~^. Filling factors for four of the minima are indicated. 

FIG. 8. The conductance characteristics G{Vsd) in a perpendicular magnetic field at T = 0.1 K, 
taken at a matched carrier density of n = 0.95 x 10^^ cm~^. Sweeps were taken at 
-6 = 0, 0.4, 0.9, 2.6 T with the filling factors indicated. The curves are offset vertically with 
zeros indicated by the dashed lines. The splitting of the zero field resonance into two peaks defines 
the gap of width A. Arrows indicate the position of inter-Landau level transitions. 

FIG. 9. The gap parameter A at matched carrier density n = 0.95 x 10^^ cm^^ as a func- 
tion of magnetic field, measured from traces similar to those shown in Fig. |^. The linear fit, 
A = 0.45^a;c — 0.19 meV, disregards measurements (open circles) where the filling factor ap- 
proaches integral or fractional (2/3) values. 
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FIG. 10. A grey-scale plot of the tunneling conductance G as a function of Vsd and mag- 
netic field B for rei = n2 = 1.48 x lO^^cm"^. Light areas correspond to regions of high conduc- 
tance. The solid lines show the splitting of the tunneling resonance, with a fitted separation of 
A = OAShuc — 0.26 meV. Dashed lines indicate inter-LL transitions at eV^d = il.OT/icJc- The 
horizontal stripes at i? = 1 T and 1.5 T show where the tunneling is completely suppressed in the 
QH regime at filling factors u = 6 and v = A, respectively. 

FIG. 11. Main figure: the gap parameter /S.y at fixed half-integral filling factors z/ as a func- 
tion of magnetic field. The solid line shows the expression = 0.45?ia;c — 0.19 meV. Inset: 
(A,, -I- 0.4 meV) 1/^/2 plotted versus B^/"^ for v > 5/2. The slope of the solid line gives o' = 2.0. 

FIG. 12. (a) The conductance characteristics GiVgd) at = 8 T and v = 1/2 at temperatures of 
T = 1.5, 1.9, 2.4, 2.8, 3.5, 4.3, 5.0, 6.0 K. (b) Arrhenius plot of G{Vsd = 0) versus 1/T at S = 8 T 
and u = 1/2. The linear fit at low temperatures gives an activation energy £'a = 0.35 it 0.01 meV. 

FIG. 13. The measured activation energy at = 1/2 (circles) and the best fit value of Eq 
(squares) , function of magnetic field. 

FIG. 14. (a) Conductance curves calculated using the double-Gaussian spectral density (Eq. |2^ 
with £^0 = 2.6 meV and a = 5.3. (b) Arrhenius plot of the calculated conductance G{Vsd = 0), 
showing an activation energy of 0.35 meV. 
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